Schrodinger cat state of a Bose-Einstein condensate in a double-well potential 
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We consider a weakly interacting coherently coupled Bose- 
Einstein condensate in a double-well potential. We show by 
means of stochastic simulations that the system could possibly 
be driven to an entangled macroscopic superposition state or 
a Schrodinger cat state by means of a continuous quantum 
measurement process. 

03.75.Fi,03.65.Bz,05.30.Jp 



I. INTRODUCTION 

Since the first observations of Bose-Einstein condensa- 
tion in dilute alkali-metal atomic gases jl]-[| the ultra- 
cold atomic gases have stimulated significant theoretical 
and experimental interest [Q . The scientific progress has 
been rapid and examples of recent experiments include 
the development of accurate detection methods the 
state preparation of topological structures 0], and the 
applications in nonlinear atom optics 0. Due to the 
macroscopic quantum coherence Bose-Einstein conden- 
sates (BECs) could possibly be also used in the future 
as a test for the foundations of quantum mechanics. One 
particularly puzzling and controversial issue has been the 
existence of macroscopic quantum superposition states in 
many-particle quantum systems. In this paper we pro- 
pose a method of creating the Schrodinger cat states of 
different atom occupation numbers in a weakly interact- 
ing BEC confined in a double-well potential. 

The existence of the superpositions of macroscopically 
distinguishable states in BECs has been addressed by 
several authors PHl^] . The superposition state may arise 
as the ground state of a coherently coupled BEC in a 
double- well potential |) 10 1. Under certain conditions 
it could be reached as a result of a unitary time evo- 
lution |ll[ . Previously, we proposed a method of creat- 
ing Schrodinger cat states in BECs by means of scatter- 
ing light from two BECs moving with opposite velocities 
||. The nonunitary evolution due to the detections of 
scattered photons drives the condensates to macroscopic 
quantum superposition states. In this paper we show 
that a continuous quantum measurement process could 
also drive a trapped coherently coupled BEC in a double- 
well potential to a Schrodinger cat state. The advantage 
of the proposed scheme is that the BEC is almost station- 
ary and trapped. Moreover, as a result of the back-action 
of quantum measurement process the superposition state 
could be reached rapidly unlike in a slow unitary evolu- 
tion, which may be very sensitive to decohcrence. 



The paper is organized as follows: We begin in Sec. H A 
by introducing the unitary system Hamiltonian. The 
scattering of li ght a nd the measu rement geometry is de- 
scribed in Sec. If B. In Sec. II C we study the dynamics 



of the open quantum system in terms of stochastic tra- 
jectories of state vectors. The results of the numerical 
simulations are presented in Sec. [Lfl|. Finally, a few con- 
cluding remarks are made in Sec. [TV|. 



II. SYSTEM DYNAMICS 

A. Unitary evolution 

We consider the evolution of a BEC in a double-well 
potential in a two-mode approximation. Macroscopic 
quantum coherence of BECs results in coherent quantum 
tunneling of atoms between the two modes representing 
'two BECs'. This is analogous to the coherent tunnel- 
ing of Cooper pairs in a Josephson junction |l^-|T8[]. To 
obtain the system Hamiltonian in the two-mode approx- 
imation for the unitary evolution of the BEC we approx- 
imate the total field operator by the two lowest quantum 
modes ip(r) ~ ipb( r )b + ip c (r)c, where ipb and ip c stand 
for the local mode solutions of the individual wells with 
small spatial overlap. The corresponding annihilation op- 
erators are denoted by b and c. The Hamiltonian in the 
two- mode approximation reads fll4j: 



«[(6t) 2 6 2 + ( c t) 2 c 2 ] 



(1) 



Here £ is the energy difference between the modes. The 
tunneling between the two wells is described by f2, which 
is proportional to the overlap of the spatial mode function 
of the opposite wells. The short-ranged two-body inter- 
action strength is obtained from k = 2-Kahfm J |V>;,(r)| 4 , 
where a and m denote the scattering length and the 
atomic mass, respectively. For simplicity, here we have 
assumed that J | V^fo C r ) 1 4 = J" |'0 c (r)| 4 . A necessary con- 
dition for the validity of the single-mode approximation 
in a harmonic trap is that the oscillation energy of the 
atoms does not dominate over the mode energy spacing 
of the trap. 

According to the Josephson effect, the atom numbers 
of the BECs determined by the Hamiltonian ([!]) may os- 
cillate even if the number of atoms in each well is initially 
equal. Due to the nonlinear self-interaction the num- 
ber oscillations also exhibit collapses and revivals. These 
have been studied numerically in Ref. p3. We may also 
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obtain a simple analytical description by solving the dy- 
namics in the rotating wave approximation in the limit 
£1 ^> Nk as described in Ref. |l9) . Here N denotes the 
total number of atoms. In particular, we may solve the 
number of atoms Nb = (b*b) in well b. We consider a 
coherent state in the both wells as an initial state. Then 
we obtain 

N b = y [1 + e^W"''- 1 ' (y/1 - (3 2 cos i] - (3 sin tp sin t])] , 

(2) 

with ij = N/3 sm(nt) cos ip — 2Q,t. Here all the operators 
on the right-hand side have been evaluated at t = 0. It 
is useful to define the real expectation values (3 and ip in 
the following way: 

(3e^ = ^c). (3) 

For a coherent state with equal atom numbers in the two 
wells we obtain the visibility (3=1. The relative phase 
between the wells is tp. For a number state there is no 
phase information and f3 = 0. For unequal atom numbers 
the maximum visibility is /3 max = 2(A r f,A r c ) 1 / 2 /7V. We see 
that the number of atoms in Eq. (Q) may oscillate in the 
case of initially equal atom numbers (3=1. The am- 
plitude of sinusoidal oscillations, representing the macro- 
scopic coherence, collapses. For instance, for ip = tt/2 
and (3=1 we may obtain the short time decay by con- 
sidering the time scales Nut <1< Vtt. Then the decay 
of the oscillations has the form exp(— NK 2 t 2 /2). This is 
the rate of the phase diffusion and it may be interpreted 
as the width of the relative phase (Aip 2 (t)) ~ Nn 2 t 2 ~ 
K 2 t 2 1 (A(p 2 (0)). Perhaps surprisingly the functional de- 
pendence of the width in this case is the same as in the 
case of two uncoupled BECs ||. 

If the phase is unknown we may obtain the ensem- 
ble average by integrating over the relative phase ip in 
Eq. @ : 

m = y {i + e McosM)-i] ^rrrpi 

x cos(2fit) J [N(3 sin(Kt)] } , (4) 

where Jo is the 0th order Bessel function. If the both 
wells have initially equal number of atoms, (3=1, the 
atom numbers do not oscillate. For unequal atom num- 
bers the oscillations collapse at the first zero of the Bessel 
function t ~ 2A/{N(3k) for (N b N c )^ 2 > 1. 

B. Quantum measurement process 

The time evolution of the system is nonunitary, when 
we include the effect of quantum measurement process. 
We consider the nondestructive measurement of the num- 
ber of atoms in the both wells by means of shining coher- 
ent light beams through the atom clouds. The scattered 



light beams are combined by a 50-50 beam splitter. We 
display the measurement setup in Fig. [I]. 




FIG. 1. The measurement setup. Two incoming light 
fields are scattered from two coherently coupled BECs. The 
two scattered photon beams are combined by a 50-50 beam 
splitter. The photons are detected from the two output chan- 
nels of the beam splitter. One of the output channels intro- 
duces only a constant phase shift and it may be ignored. 

We assume that the incoming light fields are detuned 
far from the atomic resonance. For instance, if the shape 
of the gas is flat and the light is shone through a thin 
dimension, the multiple scattering is negligible and the 
sample can be considered optically thin. A BEC atom 
scatters back to the BEC via coherent spontaneous scat- 
tering, stimulated by a large number of atoms in the 
BEC. Coherently scattered photons are emitted into a 
narrow cone in the forward direction. By spontaneous 
scattering we mean that the emission is not stimulated by 
light, although it is stimulated by atoms. The decay into 
noncondensate center-of-mass states is also stimulated by 
the Bose-Einstein statistics. However, at very low tem- 
peratures this stimulation is much weaker because most 
of the particles are in the BEC. As a first approximation 
we ignore the scattering from and to the noncondensate 
modes. Then the measurement is nondestructive in the 
sense that BEC atoms in the modes b and c scatter back 
to the same modes b and c. Because the overlap of the 
mode functions of the different wells is assumed to be 
small, the scattering between the two wells is ignored. 

The detection rate of photons on the detectors is the 
intensity of the scattered light I(r) integrated over the 
scattering directions divided by the energy of a photon 
Tick. Here k and c stand for the wave number and the 
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velocity of light, 
channel j : 



We obtain the detection rate at the 



7j 



± J d n A r 2 I j (v) = 2T{C}C j 



(5) 



The photon annihilation operator at the output channel 
j of the beam splitter is denoted by Cj . For a symmetric 
measurement geometry we obtain 



Ci 



1 



V2 



(tfb - c t c ), c 2 



1 



v/2 



(6) 



Because the total number of atoms is assumed to be con- 
served, the operator N = b^b+c'c contributes to the mea- 
surements only through a constant phase shift. There- 
fore, we may ignore the effect of the scattering channel 2 
on the dynamics. 

The scattered intensity may be written in terms of the 
positive frequency component of the scattered electric 
field E+(r) 



/(r) = 2ceo(E-(r)-E + (r)> 



(7) 



Here eo denotes the permittivity of the vacuum. 

We assume that the driving electric fields may be ap- 
proximated by plane waves E^(r) = £e e i ( k ' r-f2 *)/2. In 
the limit of large atom-light detuning A we use the first 
Born approximation and write the electric fields in the 
far radiation zone (kr ^ 1). Then the scattered field 
from the well b has the following form [f20||: 



E+(nr) 



k 2 TZe lkr 



n x (n x d) 
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dVe i(k - fcfi > r '|Vb(r')|V&. 



(8) 



Here we have defined the Rabi frequency TZ of the atomic 
dipole matrix element d by 7Z = d£/(2K). We also as- 
sumed that d • e = d. In the limit that the characteristic 
length scale I of the BECs is much larger than the in- 
verse of the wave number of the incoming light £ 3> 1/fc, 
the momentum of the scattered photon is approximately 
conserved, and we obtain in Eq. (^): 



/ 



z(k— ki) 



Vb(r')| 2 -S(kh-k) 



In this simple case the scattering rate T may be easily 
evaluated: 



r = 



3 7 ^ 2 
8ttA 2 ' 



(9) 



Here 7 = d 2 k 3 / (6irheo) denotes the optical linewidth of 
the atom. 



C. Stochastic Schrodinger equation 

The dissipation of energy from the quantum system 
of macroscopic light fields and the BEC in a double-well 
potential is described by the coupling to a zero temper- 
ature reservoir of vacuum modes, resulting in a spon- 
taneous emission linewidth for the atoms. The dynam- 
ics of the continuous quantum measurement process may 
be unraveled into stochastic trajectories of state vec- 
tors [pT| p3|. The procedure consists of the evolution 
of the system with a non-Hermitian Hamiltonian H e g, 
and randomly decided quantum 'jumps'. In our case the 
quantum jumps correspond to the detections of sponta- 
neously emitted photons. The system evolution is thus 
conditioned on the outcome of a measurement. The non- 
Hermitian Hamiltonian has the following form: 



H eS = H S - ihTCld , 



(10) 



where the unitary system Hamiltonian H$ is determined 
by Eq. ®. 

Equation ( |io| ) corresponds to the modification of the 
state of the system associated with a zero detection re- 
sult for scattered photons. Because the output is being 
continuously monitored, we gain information about the 
system even if no photons have been detected. 

The Hamiltonian H c s determines the evolution of the 
state vector ip sys (t). If the wave function ip sys (t) is nor- 
malized, the probability that a photon from the output 
channel 1 of the beamsplitter is detected during the time 
interval [t, t + 5t] is 



P(t) = 2T(ij sys (t) I Cld I Vsys(t)) St . 



(11) 



We implement the simulation algorithm as follows: At 
the time to we generate a quasi-random number e which is 
uniformly distributed between and 1. We assume that 
the state vector ip S ys(to) at the time to is normalized. 
Then we evolve the state vector by the non-Hermitian 
Hamiltonian H c g iteratively for finite time steps At. At 
each time step n we compare e to the reduced norm of 
the wave function, until 

(ip sys (t + nAt) I 4> sys {t + nAt)) < e , 

when the detection of a photon occurs. If the photon has 
been observed during the time step t — > t + At we take 
the new wave function at t + At to be 



(12) 



\4> sys (t + At)) = V2rCi IVsysW), 

which is then normalized. 

III. NUMERICAL RESULTS 



We simulate the effect of the system Hamiltonian H$ 
and the quantum measurement process of scattered pho- 
tons by means of the stochastic Schrodinger equation. 
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For simplicity, we set the total number of atoms to be 
reasonably small N = 200. We start from a slightly 
asymmetric initial state with the two modes in number 
states N b = 102 and N c = 98. We choose T/Q = 5x 10~ 6 . 

After just a few detected photons we observe the emer- 
gence of two well-separated amplitude maxima in the oc- 
cupation number of atoms in one of the two wells. These 
correspond to a macroscopic number state superposition 
or a Schrodinger cat state. Because the total number of 
atoms is assumed to be conserved, the atom numbers in 
the two wells are entangled and we have a Bell-type of su- 
perposition state. In Fig. ||we display the absolute value 
of the wave function \tpb\ in mode b in number state basis 
in a single run at two different times. In this case the 
nonlinearity vanishes n — and £/f2 = 0.1. We clearly 
recognize the two distinct peaks in the number distribu- 
tion. For instance, the peaks in Fig. || (b) are centered 
at Ni, ~ 10 and Nb — 190 corresponding to maxima at 
iV c ~ 190 and N c ~ 10, respectively. In Fig. || we show 
the absolute value of the wave function for a different run 
with Nk/Q = 0.2 and £/fi = 0.001. 

We also describe the state of the BEC in terms of the 
quasiprobability Q distribution. For the number state 
distribution of atoms \ipb) = Cn\n) in mode b we ob- 
tain pi: 



Q(a) 



\(a\ipb 



-M 2 



N 

E 



arc: 



(13) 



The Q function represents the phase-space distribution. 
The amplitude and phase quadratures are denoted by X 
and Y. In polar coordinates the radius in the xy plane 

1/2 

is equal to N b and the polar angle is the relative phase 
between the atoms in the two wells. In Fig. ^ we show the 
Q function distribution of the number state superposition 
displayed in Fig. |2|(a). 

It is interesting to emphasize that the measurement of 
the number of atoms in only one of the wells affects the 
system dynamics quite differently. In that case the num- 
ber state distribution remains well localized and approxi- 
mately approaches a coherent state 17 . Even though we 



start from a number state with no phase information, the 
detections of spontaneously scattered photons establish 
a macroscopic coherence or the off-diagonal long-range 
order (ODLRO) between the atoms in the two separate 
wells. This is similar to establishing the coherence be- 
tween two BECs as a result of the counting of atoms 
|p5| , p6| . However, in the present case the continuous mea- 
surement process drives the system to a Schrodinger cat 
state and the ODLRO remains small. We may describe 
the visibility of the macroscopic coherence between the 
two wells by the real parameter defined in Eq. (||) . We 
show the relative visibility (3 r = /3//3 max and the number 
of atoms in well b as a function of the number of measure- 
ments for k = in Fig. || and for Nn/fl = 0.2 in Fig. 
Due to the emergence of the superposition state the vis- 
ibility remains below one. The measurement process of 



the scattered photons significantly complicates the dy- 
namics of the number of atoms predicted by the unitary 
time evolution of Eq. (h]) . 





FIG. 2. The entangled Schrodinger cat states of atoms 
with different atom numbers. We show the absolute value of 
the wave function |i/>{,| in the number state basis for the atoms 
in well b during one realization of stochastic measurement 
process of spontaneously scattered photons after (a) 750 and 
after (b) 1700 detections. The two maxima correspond to the 
superposition states. The nonlinearity k = and the total 
number of atoms N = 200. 




FIG. 3. The entangled Schrodinger cat states of atoms 
with different atom numbers. We show the absolute value 
of the wave function during one realization of stochastic 
trajectory after (a) 100 and after (b) 1250 detections. The 
nonlinearity Nn — 0.2. 




FIG. 4. The Q function of the Schrodinger cat state. We 
show the Q function of the quantum state displayed in Fig. ^ 
(a). The two peaks are located at different value of the ra- 

1/2 

dius N b representing the different maxima of the occupation 
numbers. 
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FIG. 5. The (a) relative visibility of the interference j3 r 
and the (b) number of atoms Nt in one of the wells as a 
function of the number of detected photons. Due to the 
Schrodinger cat state the visibility never reaches one. Here 
the nonlinearity k = 0. 
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FIG. 6. The (a) relative visibility of the interference j3 r 
and the (b) number of atoms Nt in one of the wells as a 
function of the number of detected photons. The nonlinearity 
Nn/fl = 0.2. 



IV. FINAL REMARKS 

We studied the generation of the macroscopic super- 
position states or the Schrodinger cat states of a BEC in 
a double-well potential. The Schrodinger cat state was 
shown to emerge as a result of the continuous quantum 
measurement process of scattered photons. The partic- 
ular detection geometry increses the fluctuations of the 
relative atom number between the two wells. Therefore 
the superposition states are more stable in the detection 
process. The proposed setup is an open quantum system 
and the creation of the Schrodinger cat state in this case 
is not based on reaching the ground state of a BEC in a 
double- well potential [p|JTo| . The advantage over previ- 
ously proposed open systems schemes |1 is that the BEC 
is stably trapped and the superposition state for a small 
BEC could be created by scattering only a few photons. 

In the present discussion we ignored the effect of deco- 
herence f27j. The interaction of the BEC with its envi- 
ronment results in the decoherence of the superposition 
states. We can identify several sources of decoherence. 
Decoherence by amplitude damping or by phase damping 
has been estimated in Ref. p8[ . The inelastic two-body 
and three-body collisions between the condensate atoms 
and the noncondensate atoms change the number of con- 
densate atoms and introduce amplitude damping. The 



phase damping corresponds, e.g., to elastic collisions be- 
tween the condensate and noncondensate atoms in which 
case the number of BEC atoms is conserved. If the num- 
ber of atoms in a BEC is not large, the scattering between 
the condensate and noncondensate atom fractions may 
not be negligible. This also introduces amplitude deco- 
herence. Additional sources of decoherence may be, e.g., 
the imperfect detection of the scattered photons and the 
fluctuations of the magnetic trap. In Ref. Q it was pro- 
posed that the decoherence rate of a BEC could be dra- 
matically reduced by symmetrization of the environment 
and by changing the geometry of the trapping potential 
to reduce the size of the thermal cloud. Moreover, the 
continuous measurement process increases the informa- 
tion about the system and therefore it could also reduce 
the decoherence rate. 
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